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SECTION  1  -  INTRODUCTION 


Previous  studies  (References  1  and  2)  have  Indicated  that  a  quadrilateralized 
spherical  cube  dissection  of  the  Eardi  will  provide  a  uniquely  simple  and  useful 
set  of  coordinates  in  \^ich  to  store  remote  sensor  data  on  a  global  scale.  This 
representation  divides  the  globe  into  six  square  grids,  allowing  an  extremely 
simple  and  elegant  addressing  scheme. 

A  mapping  fimction  (Reference  1)  has  been  developed  which  permits  an  equiva¬ 
lent  (equal  area)  mapping  of  a  sphere  onto  the  faces  of  a  cube.  This  mapping 
provides  an  ideal  match  to  constant  resolution  sensors  and  allows  the  simplifi¬ 
cation  of  map  factor  calculations  for  meteorological  applications.  Figure  1-1 
illustrates  the  cubic  dissection  of  a  sphere,  and  Figure  1-2  shows  a  map  of  the 
world  as  projected  using  the  equivalent  mapping  function. 

In  studying  the  feasibility  of  Implementing  this  quadrilateralized  spherical 
cube  Earth  data  base  system,  a  number  of  areas  were  identified  which  were 
deemed  worthy  of  further  study.  In  the  course  of  these  studies,  additional 
details  and  problems  were  brought  to  light.  The  initial  objectives  of  this 
study  were 

•  Improvement  of  the  formulation  of  the  data  base  mapping  functions 

•  Study  of  methods  of  representing  satellite  motions  in  terms  of  the 
data  base  coordinate  system 

•  Study  of  data  base  Input/Output  optimization  on  the  CYBER/175 
computer 

•  Production  of  maps  illustrating  various  geographic  orientations  of 
the  data  base  cube 
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Figure  1-1.  Cubic  Dissection  of  a  Sphere  by  an 
Inscribed  Cube 


1-2 


Additional  studies  required  were 


•  Study  of  algorithms  required  to  implement  a  revised  data  base  fill¬ 
ing  scheme 

•  Improvement  in  the  accuracy  of  the  data  base  mapping  functions 
and  the  matching  of  direct  and  Inverse  functions 

The  general  purpose  in  these  studies  has  been  twofold:  (1)  to  optimize  the 
existing  data  base  scheme  and  (2)  to  assist  in  the  timely  development  of  a  pre¬ 
liminary  implementation  by  in:^)rovlng  and  clarifying  previously  developed 
formulae  and  algorithms.  Most  sections  of  this  report  deal  with  the  mainte¬ 
nance  of  current  data  in  the  data  base.  The  portion  of  the  system  which  will 
fulfill  these  functions  has  been  designated  the  preprocessor  program. 

This  report  presents  several  new  and  refined  algorithms  and  mappings  aimed 
at  improving  the  speed  and  accuracy  of  an  operational  preprocessor  program. 
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SECTION  2  -  BACKGROUND 


2.1  PRIOR  RESULTS 

A  proposed  set  of  algorithms  for  implementing  the  quadrilateralized  spherical 
cube  data  base  was  presented  in  another  report  (Reference  1).  That  report  ad¬ 
dressed  numerous  problems:  mapping  points  from  the  Earth's  surface  into  the 
data  base,  determining  serial  addresses  for  data  base  locations,  fast-filling 
and  data  retrieval,  and  determining  the  ground  locations  of  scanner  pixels.  The 
results  presented  herein  are  derived  from  further  study  of  these  problems. 

For  background  information.  Sections  2.2  and  2.3  present  the  previously  de¬ 
rived  transformation  and  method  of  fast-filling. 

2.2  PREVIOUSLY  DEVELOPED  TRANSFORMATION 

Section  2  and  Appendix  A  of  Reference  1  present  the  mathematical  development 
of  the  mapping  function,  f  ,  and  its  inverse,  f*  ,  These  functions  relate  the 
data  base  coordinates  (x  ,  y)  to  the  rectangular  coordinates  ,  rj)  of  the 
point  on  any  cube  face  resulting  from  the  central  projection  of  a  datum  point 
on  the  sphere.  For  convenience,  the  mapping  functions,  f  and  f*  ,  are  re¬ 
peated  here.  The  coordinates  ,  tj)  are  related  to  the  data  base  coordinates 
(X  ,  y)  via 


i  =  f(x,  y) 
T?  =  f(y.  X) 


(2-1) 
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where  the  function  f  ,  called  the  direct  hinctlon,  is  given  by 


where 


f(x,  y)  =  yx  + 


(1-y) 


o 


3 

X 


(2-2) 


6  »  ~  [-(M  +  2y)  +  -  44y  +  4y^  +  16  v^y^)] 

4r'* 

o 

»  0. 79048  64491  208 

iji  =  —ij  (3  -  2y  -  fl  -  2r^  sj  =  -1. 2264  41487  984 

2r^  '  ° 

o 


y  =  v/irTe 


r  =l/v/3 
o 


(2-3a) 


(2-3b) 


(2-3c) 

(2-3d) 


(2-3e) 
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and  the  coefficients 


'ij 


and  d^  are 


^00  = -2.7217 

05366 

1814 

d  =  1.4833 
0 

12929 

4187 

=  -5.5842 

16830 

5430 

d^  =  1. 1199 

72606 

9742 

01 =  2. 1711 

17480 

9423 

d^  =  6.0515 

38216 

1464 

:20  = -3.4578 

62747 

3390 

=  -6.4160 

15152 

6783 

02  =  1.9736 

26575 

8872 

The  inverse  function,  f*  ,  relates  (x  ,  y)  to  (4  »  i?)  ;  that  is. 


x  =  f*  (4,  TJ) 

y  =  f*  (t).  i) 


(2-4) 


and  is  given  by 


V)  = 


2  ^ 
r 
o 


2i^2i 


(2-5) 
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where 


y*  =  i/y 
H*  =  i/n 


yj  =  i/(y  + 

Mj  =  l/(M  +  2r^6) 
6*  =  (M*  -  H*)/2r‘^ 


(2-6a) 

(2-6b) 

(2-6c) 

{2>6d) 

{2-6e). 

(2-6f) 

(2-6g) 


and  the  coefficients  c*  and  d*  are 

Cqq*  =  3.973  89249 
Cio*  =  6.591  19476 
Coi*  =  -25.368  92536 
C2o*  = -73.064  97000 
c^^*  =  77.381  61133 
Cq2*  =  21.685  89623 


dp*  “  1.811  28250 
d^*"  37.635  47857 
dg*"  63.000  23655 


The  use  of  f*  in  an  algorithm  for  obtaining  the  data  base  coordinates  (x  ,  y) 
of  any  datum  point  on  the  Earth's  surface  is  described  in  Section  3. 3.  The 
algorithm  presented  in  Section  3. 3  is  the  result  of  the  present  study  to  obtain 
a  method  of  implementing  f  *  in  a  con^utationally  efficient  manner. 
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2.  3  PREVIOUSLY  DEVELOPED  FAST- FILLING  ALGORITHM 


Sections  4  and  5  of  Reference  1  present  a  "fast-filling' ’  method  for  updating  the 
data  base.  The  method  was  developed  to  obtain,  as  directly  as  possible,  the 
data  base  location  corresponding  to  the  location  of  a  given  scanner  data  element 
and  to  store  the  value  associated  with  that  element  in  the  calculated  data  base 
location.  Starting  with  a  given  satellite  location  and  scanner  angle,  a  series 
of  five  steps  are  necessary  to  calculate  the  corresponding  data  base  location. 

1.  Determine  the  three-dimensional  Cartesian  coordinates  of  the  point 
observed  by  the  scannei'. 

2.  Determine  the  face  onto  which  the  point  is  to  be  mapped 

3.  Project  from  the  Cartesian  (a,  b,  c)  system  into  an  intermediate 
system  (4  ,  7?)  as  shown  in  Figure  2-1 

4.  Transform  from  ,  17)  into  two-dimensional  Cartesian  coordi¬ 
nates  (X  ,  y)  on  the  face  of  the  data  base  cube 

5.  Convert  from  face  number  and  (x  ,  y)  to  a  serial  location  in  the 
data  base 

This  transformation  need  not  be  performed  for  every  point,  since  both  the  pro¬ 
jection  onto  the  Earth's  surface  and  the  (( ,  17)  to  (x  ,  y)  transformation  may 
be  linearly  approximated  over  a  relatively  long  range.  Thus,  the  complete  com 
putation  is  required  only  for  widely  spaced  benchmark  points,  with  intermediate 
points  being  determined  by  interpolation. 

Input  and  output  are  handled,  as  described  in  Section  5  of  Reference  1,  by 
blocking  both  the  input  data  and  the  data  base  into  rectangular  arrays.  Methods 
for  handling  problems  such  as  queuing  data  base  I/O,  splitting  ii4)ut  blocks 
across  face  boimdaries,  and  locating  a  vertex  within  an  input  block  were  also 
discussed  in  Reference  1. 
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In  Reference  3,  this  scheme  was  compared  with  an  alternative  scheme.  The 
alternative  scheme  follows  the  above  method  through  step  3;  then,  instead  of 
storing  from  the  input  block  directly  into  the  data  base,  the  following  steps  are 
performed: 

•  For  each  input  data  block,  determine  which  data  base  blocks  will  be 
completely  covered  by  the  input  data 

•  Determine,  in  the  ,  77)  system,  a  set  of  benchmark  points  cor¬ 
responding  to  the  comers  of  the  data  base  blocks  to  be  filled 

•  Use  the  above  benchmarks  to  interpolate  the  ,  ?))  coordinates 
corresponding  to  each  (x,  y)  data  base  location  and  store  the  closest 
point  in  the  input  data  block  into  the  corresponding  serial  location 

This  method  offers  the  advantage  of  providing  a  completely  filled  data  base 
block  and  eliminates  data  base  read  accesses  during  the  update  process.  After 
a  thorough  examination  of  this  method,  CSC  concluded  that  the  advantages  out¬ 
weigh  the  slight  addHional  computational  cost  incurred  in  finding  the  appropriate 
data  base  blocks.  Consequently,  some  of  the  work  of  this  report  was  oriented 
toward  this  approach. 
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SECTION  3  -  NEW  PROJECTION  AND 
TRANSFORMATION  METHODS 


3. 1  GROUND  TRACE  EQUATIONS  IN  KEPLERIAN  ELEMENTS 

This  section  presents  general  relations  that  yield  the  geocentric  coordinates  of 
the  subsatellite  point  and  of  any  point  on  the  Earth's  (spherical)  surface  lying 
on  an  arc  scanned  by  an  onboard  satellite  sensor.  These  geocentric  coordinates 
are  to  be  computed  as  an  intermediate  step  in  algorithms  for  computing  data 
base  coordinates  (Sections  3.2  and  3.3). 

The  following  subsections  distinguish  between  two  cases:  (a)  the  computation 
of  the  coordinates  of  the  subsatellite  point,  and  (b)  the  computation  of  the  coor¬ 
dinates  of  any  point  on  the  scan  arc.  Although  the  first  case  could  be  treated 
as  a  special  case  of  the  second.  Its  development  is  separated  in  order  to  facil¬ 
itate  computer  implementation.  The  following  subsections  sununarize  the 
results,  which  are  developed  in  detail  in  Appendix  A. 

3. 1. 1  Subsatellite  Point 

It  is  assumed  that  the  satellite  ephemeris  is  available  in  the  form  of  Keplerian 
orbital  elements.  Expressions  can  be  obtained  in  terms  of  these  elements  for 
the  Earth-centered  direction  cosines  (1  ,  m  ,  n)  of  the  subsatellite  point  in 
the  equatorial  x,  y,  z  coordinate  system  shown  in  Figure  3-1. 

Since  the  subsatellite  point  lies  on  the  geocentric  position  vector  of  the  satellite 
(see  Figure  3-1),  1  ,  m  ,  and  n  will  be  identical  with  the  corresponding  direc¬ 
tion  cosines  for  the  satellite's  orbital  position.  The  desired  expressions  for  1  , 
m  ,  and  n  are  developed  in  Appendix  A  to  be 

1  =  +COS  n  cos  (to  +  f)  -  sin  O  sin  (to  +  f)  cos  i 
m  =  +sin  n  cos  (to  +  f)  +  cos  O  sin  (to  +  f)  cos  i  (3-1) 

n  =  +sin  (to  +  f)  sin  i 
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Figure  3-1.  Keplerian  Orbital  Elements  and  Geographic  Coordinates 
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The  angular  arguments  are  Kepleiian  elements  and  are  referenced  to  the  true 
equator  and  prime  (Greenwich)  meridian  of  date  coordinate  system.  The  di¬ 
rection  cosines  (1  ,  m  ,  n)  in  Equation  (3-1)  are  obtained  for  the  subsatellite 
point  corresponding  to  any  point  in  the  orbit  using  the  true  anomaly,  f  ,  which 
reflects  the  motion  of  the  satellite  in  its  orbit,  and  the  time-dependent  node,  O  , 
which  reflects  the  rotation  of  the  Earth.  The  time  dependence  of  the  node  is 
given  by 


0  =  0  -  s  (3-2) 

o 

where  O^  is  the  initial  (constant)  value  of  O  and  s  is  the  elapsed  Greenwich 
sidereal  time  since  the  instant  that  0  =  0^.  For  high  precision,  (* ,  a; ,  i  , 
and  s  should  be  corrected  for  the  Earth's  nutation,  precession,  and  polar 
motion. 

3.1.2  Nonzero  Scan  Angle 

It  is  assumed  that  the  satellite  collects  data  through  an  onboard  sensor  which 

scans  the  surface  of  the  Earth  through  an  arc  of  fixed  size  in  such  a  way  that  its 

sweep  carries  it  throi^h  the  subsatellite  point  in  a  direction  at  right  angles  to 

the  direction  of  the  satellite's  orbital  motion.  The  trace  of  the  scan  arc  on  the 

surface  of  the  (spherical)  Earth  is  a  segment  of  a  great  circle  subtending  a 

fixed  central  angle,  ,  as  illustrated  in  Figure  3-2.  The  introduction  of  the 

central  angle,  ip  ,  though  generally  applicable  in  defining  the  groimd  location 
o 

of  a  scan  element  (datum  point),  is  particularly  suitable  when  scan  elements 
are  spaced  equally  along  the  great  circle  arc. 

Expressions  are  presented  below  for  obtaining  the  geocentric  longitude,  0'  , 
and  colatitude,  0'  ,  of  any  pair  of  points  (datum  points)  on  the  scan  arc  which 
are  located  symmetrically  with  lespect  to  the  subsatellite  point.  Details  of 
their  development  are  given  in  Appendix  A.  These  expressions  are  then 
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z 


SUBSATELLITE  POINT 


SCAN  ARC  TRACE 


♦  V 


ORBITAL  GROUND 
TRACE 


Figure  3-2.  The  Scan  Arc  as  a  Central  An^e 
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employed  in  Sections  3. 2  and  3.3  in  describing  the  two  alternative  algorithms 
for  obtaining  the  data  base  coordinates  (x  ,  y)  of  the  datum  points. 

First,  the  sines  and  cosines  of  the  geocentric  longitude,  6  ,  and  colatitude,  0  , 
of  the  subsatellite  point  are  required  in  terms  of  the  direction  cosines  available 
from  Equation  (3-1),  i.e. , 


sin  0  cos  6  =  1  (3-3a) 

sin  0  sin  0  =  m 

cos  0  =  n 

Actually,  Equation  (3-3c)  yields  both  cos  0  and  sin  0  ,  since  0  ^  0  ^  180  de¬ 
grees.  Consequently,  sin  0  and  cos  0  can  be  determined  from  Equations  (3-3a) 
and  (3-3b). 

Returning  to  the  scan  arc,  it  is  convenient  to  view  (()/2  as  a  general  variable 

t 

describing  the  positions  of  any  pair  of  points  along  tlie  arc.  The  following  ex¬ 
pressions  yield  two  unique  values  for  both  6*  and  0'  ,  depending  on  the  alge¬ 
braic  sign  of  )/)/2  .  Thus, 

cos  0'  =  cos  0  cos  il>/2  -  cos  i  sin  il>/2  (3-4) 

where  cos  0  and  cos  i  are  already  available.  Equation  (3-4)  gives  0' 
uniquely,  since  0^0'^  180  degrees,  from  which  sin  0'  is  then  computed 
and  used  in 


(3-3b) 

(3-3c) 


sin  (6'  -  6)  = 


cos  (to  +  f)  sin  i  sin  il)/2 
sin  0  sin  0' 


(3-5) 
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(3-6) 


cos  (6*  -  6)  = 


cos  0/2  -  cos  0  cos  0* 
sin  0  sin  0* 


to  obtain  6'  ,  where  all  of  the  terms  to  the  right  of  the  equalities  have  pre¬ 
viously  been  determined.  • 

Caution  is  necessary  when  applying  Equations  (3-5)  and  (3-6)  near  the  poles, 
since  these  relations  are  singular  for  subsatellite  points  or  datum  points  near 
the  poles.  This  is  expected,  since  the  geocentric  longitudes  6  and  6'  are 
undefined  at  the  poles.  No  problem  need  arise  in  computations,  however,  since 
the  satellites  are  not  expected  to  have  exactly  90-degree  inclinations  and  since, 
for  the  case  of  an  equatorial  orientation  of  the  data  base  cube,  it  is  a  simple 
matter  to  regard  the  cube  face  coordinates  of  such  datum  points  as  zero  within 
an  aibitrary,  but  nonzero,  distance  from  the  pole. 

It  is  emphasized  that  Equations  (3-4),  (3-5),  and  (3-6)  yield  two  unique  values 
for  both  6’  and  0'  ,  depending  on  the  choice  of  ±(|)/2  ;  therefore,  both  halves 
of  the  scan  arc  are  treated  efficiently  using  the  same  auxiliary  quantities  (sin  i  , 
sin  0  ,  etc.). 

The  expressions  just  presented  will  be  used  in  Section  3. 2  in  an  algorithm  em¬ 
ploying  a  newly  developed  transformation  to  obtain  data  base  coordinates  of 
datum  points.  These  results  are  also  used  in  Section  3. 3  in  an  algorithm  de¬ 
signed  to  employ  the  older  transformation  of  Section  2. 2  to  compute  the  data 
base  coordinates  of  datum  points  in  a  computationally  efficient  manner. 

3. 2  POLAR  FORMULATION  OF  THE  AREA-PRESERVING  TRANSFORMATION 

Section  2  of  Reference  1  presents  an  area-preserving  transformation  from  the 
spherical  square  to  the  inscribed  plane  square.  That  transformation  utilizes  a 
central  projection  onto  the  plane  square  and  a  transformation  from  the  projected 
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coordinates  (^,  rf)  to  the  data  base  coordinates  (x,  y)  .  This  transformation, 
designated  by  f*,  is  given  in  Equation  (2-4)  as 

X  =  f*(4,  n) 
y  =  Ht),  0 


and  is  based  on  the  Jacobian  determinant.  It  was  noted  in  Reference  1  that  the 
solutions  obtained  to  differential  equations  based  on  this  determinant  are  not 
fully  constrained  by  requiring  the  interior  and  edges  of  the  spherical  square  to 
map  onto  the  plane  square.  Additionally,  it  was  found  that,  despite  adding  fur¬ 
ther  constraints,  no  analytic  solution  could  be  found  for  f*.  An  approximate 
solution  was  arrived  at  using  a  numerical  steepest  descent  method.  Under  the 
current  contract,  an  alternative  solution  was  studied;  the  results  are  presented 
below. 

Two  sets  of  polar  coordinate  systems  are  considered  (Figures  3-3  and  3-4). 
Figure  3-3  shows  a  spherical  square  of  radius  R,  produced  by  the  central  pro¬ 
jection  of  an  inscribed  square  upon  the  sphere.  The  polar  coordinates  (6,0) 
are  defined  by  an  axis  of  the  sphere  passing  through  the  center  of  the  spherical 
square  and  by  a  great  circle  arc  passing  through  the  axis  and  bisecting  one 
edge  of  the  square.  The  azimuthal  angle,  6  ,  is  measured  from  this  "meridian" 
and  the  central  angle,  0  ,  is  measured  from  the  axis.  For  the  case  where  the 
coordinate  axis  and  meridian  coincide  with  those  of  the  sphere,  8  equals  the 
longitude  and  0  is  the  colatitude. 

Figure  3-4  shows  a  plane  square  of  side  2r^.  The  axis  in  this  system  is  nor¬ 
mal  to  the  plane,  passes  through  the  center  of  the  plane,  and  has  its  origin  at  a 
distance  r^  below  the  plane.  Coordinates  (fx,  u)  are  defined  in  the  same  sense 
as  (6,  0)  . 
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Figure  3-3.  Polar  Coordinate  System  on  Spherical  Square 


Figure  3-4.  Polar  Coordinate  System  on  an  Inscribed  Plane 
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If  the  axes  of  the  two  systems  are  aligned  and  the  plane  square  is  inscribed 
within  the  spherical  square,  it  can  be  shown  that 


=  R/yi”  (3-7) 

The  area  of  the  sphere  is 

A  =  4irR^  (3-8) 

s 

and  the  surface  area  of  the  inscribed  cube  is 

A  =  24r^  (3-9) 

c  o 


The  ratio  of  the  areas  of  the  spherical  and  plane  squares  is  given  by 


X 


(3-10) 


Substituting  Equation  (3-7)  into  Eqiiation  (3-10)  yields 

X  =  ir/2  (3-11) 

To  further  constrain  the  mapping  of  the  spherical  square  upon  the  plane,  con¬ 
sider  the  condition  /i  =  f(6)  ,  where  n  is  Independent  of  0  .  This  condition 
implies  the  useful  property  that  any  great  circle  on  the  sphere  passing  through 
the  center  of  the  spherical  square  will  map  as  a  straight  line  passing  through 
the  center  of  the  plane  square.  This  function  must  be  area  preserving  in  the 
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sense  that  sectors  of  equal  area  on  the  sphere  must  map  into  equivalent  areas 
on  the  plane.  If  the  function  f  exists^  then  it  Is  possible  to  relate  the  areas  of 
differential  elements  d6  and  d^  of  the  spherical  and  plane  squares,  as  shown 
below. 

To  obtain  the  f\mction  ^  =  f(6)  ,  consider  the  right  spherical  triangle  ABC  in 
Figure  3-5,  and  let  one  side,  B,  be  the  meridian  and  A,  the  azimuth  angle,  6  . 
For  a  unit  sphere,  R  =  1  and  b  =  tr/4  . 

From  the  general  relation  for  the  area  of  a  spherical  triangle 

S  =  R(A  +  B  +  C  -  ff)  (3-12) 

and  the  area  of  the  right  spherical  triangle  is 

S  =  0  +  /3  -  ir/2  (c  E  tr/2)  (3-13) 

s 


By  Napiers'  identity,  B  =  cos  ^(sin  A  cos  b)  and 

S  =  0  +  cos"^(sin  0  cos  ff/4)  -  n/2  (3-14) 

s 

This  triangle  may  represent  any  segment  of  the  spherical  square  in  the  range 
-ir/4  S  0  S  ff/4  . 

The  area  of  the  corresponding  right  triangle  on  the  plane  (Figure  3-6)  is 


S 

P 


2 


(3-15) 
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Figure  3-5.  Right  Spherical  Triangle  ABC 


Figure  3-6.  Plane  Right  Triangle  Corresponding  to  ABC 
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or,  since  a  =  r^tan  n  , 


(3-16) 


For  n  =  £(6)  to  preserve  area,  the  condition  must  hold  that,  for  aiqr  0  , 


S  /S  =X 
s  p 


(3-17) 


Thus, 


S  =  ir/2  S 
s  p 


and 


In  the  case  of  a  unit  sphere,  r  =  lA/3',  and 

o 


(3-18) 


(3-19) 


(3-20) 


H  =  tan’^  ^0  -  cos^  |sin  ®  4  | |  (3-21) 

It  can  be  shown  that  n  =  0  at  0=0  and  n  =  ±ff/4  at  0  =  ±ff/4  .  The  func¬ 
tion  f  is  defined  over  the  quadrant  -ff/4  S  <0,  ir/4  and,  by  reflection, 
over  the  entire  spherical  and  plane  squares. 
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To  relate  u  and  0  ,  expressions  for  the  differential  area  elements,  0d6  and 
udn  ,  must  be  obtained.  Figure  3-7  illustrates  a  spherical  zone  of  one  base, 
having  area 


S  =27rRh 
s 

Substituting  h  =  R  -  R  cos  0  and  taking  R  =  i  yields 

Sg  =  2ir(l  -  cos  0) 


(3-22) 


(3-23) 


Since 


■r- 

J  o 


S  =  1  1  -  cos  0  d0 

s 


dS  =  (1  -  cos  0)  d6 
8 


(3-24) 


The  area  of  the  corresponding  sector  on  the  plane  is  given  by 


dS  =  r 
P 


(3-25) 


In  order  to  apply  the  area-preserving  condition,  Equation  (3-17)  is  differentiated 
to  yield 


dS  =  jr/2  dS 
s  p 


(3-26) 


Substituting  Equations  (3-24)  and  (3-25)  into  Equation  (3-26)  gives 


r  du  =  2/ir  (1  -  cos  0)  de 


(3-27) 
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Figure  3-7.  Spherical  Zone  of  One  Base 


Since  V  ,  Equation  (3-27)  can  be  rewritten  as 


2 

tan  p  = 


^(l-cos  0)|i 


(3-28) 


For  fixed  values  of  6  and  u  ,  Equation  (3-28)  becomes 


tan  =  C(1  -  cos  0) 


(3-29) 


where 


C 


ffr 


2_  ^ 
2  du 


0,u 


(3-30) 


Clearly,  C  can  be  evaluated  from  boundary  coiKUtions: 


1.  When  0  =  0,  p  =  0 


points  at  the  center  map  onto  themselves 


2.  When  0  =  0niax 
^  ~  ^max 


points  on  the  boundary'  of  the  spherical 
square  must  map  onto  points  on  the 
boundary  of  the  plane  square 
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To  develop  these  boundary  conditions.  It  is  necessary  to  obtain  0  and 

max 

p  as  functions  of  6  and  u  ,  respectively.  Figures  3-8  and  3-9  illustrate 
max 

these  angles.  It  should  be  noted  that  to  obtain  the  boundary  of  the  spherical 
squai'e,  a  central  projection  to  the  inscribed  plane  is  employed.  This  projection 
preserves  both  9  and  0  and  reduces  the  problem  to  the  case  of  the  plane 
square. 

The  loci  of  points  describing  these  edges  can  be  expressed  as 

r  =  r  sec  0  (3-31) 

8  o 


and 


r  =  r  sec  u 
P  o 


(3-32) 


and  the  central  angles  can  be  expressed  as 


0  =  tan”^  sec  0 

max 


(3-33) 


and 


P  =  tan  ^  sec  u 
max 


(3-34) 


Substituting  0  and  P  in  Equation  (3-27) 

in3X  1T13X 


tan  p  =  C(1  -  cos  0  ) 

max  max 


(3-35) 
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POINT  DEFINED  BY  (6.0} 


Figure  3-8.  Relationship  of  to  6 


Figure  3-9.  Relationship  of  to  fj. 
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Therefore 


C  = 


tan^  V  2-1 

max  _  tan  (tan  sec  u) 


2 

sec  u 


1  -  cos  ^  1  ggj,  0j  j  _  1  ggg 


(3-36) 


Substituting  C  into  Equation  (3-35)  yields 


tan^  U  = 


2 

sec  u 


1  -  cos  tan  ^  sec  6 


(1  -  cos  0) 


(3-37) 


and 


V  =  tan 


-1  /  sec^  u(l  -  cos  0) 

1  -  cos  tan  sec  8 


(3-38) 


This  satisfies  the  boundary  condition  0  =  0,  U  =  0  ,  and  additionally  has  the 

Drooertv  that  u  =  0  at  0  =  0  and  0  =  u  ,  where  0  =  iff/4  ,  u  =  iff/4  • 

f  f  j  max  max 

This  relation  holds  from  0=0  to  0  =  ff/2  ,  and  because  0  is  never 

ITictX 

greater  than  tan  ^\/2  =  0.95  ,  it  is  correct  for  the  entire  spherical  square. 
Thus,  Equations  (3-38)  and 


(3-39) 


represent  an  area-preserving  transformation  from  the  spherical  square  to  the 
inscribed  plane  square  for  the  quadrant. 

To  treat  the  remaining  quadrants,  each  is  successively  rotated  into  the  quadrant 
of  definition  by  subtracting  multiples  of  ff/2  ,  as  shown  in  Table  3-1. 
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Tab  e  3-1.  Rotations  of  6  and  u’  for  a  Given  0’^ 


Range 

0 

u' 

0'  si 

4  4 

0’ 

u 

VI 

V 

6*  -  — 

2 

u^iL 

2 

V 

<D 

V 

0’  -If 

u  +  IT 

5!L<  e- 

4  4 

6' -is 

2 

3ir 

u  +  — 

2 

a 

The  quantities  0'  and  u*  are  input  and  output, 
respectively;  and  0  and  u  are  as  used  in  the 
computation  of  Equations  (3-38)  and  (3-39). 

3. 3  REFINEMENT  OF  THE  CARTESIAN  TRANSFORMATION 

Section  4. 2  of  Reference  1  outlines  a  method  for  obtaining  the  data  base  coor¬ 
dinates  (X  ,  y)  of  any  datum  point  from  Kepleiian  orbital  elements,  time,  a 
scan  angle,  and  parameters  describing  the  attitude  of  the  satellite  using  the 
inverse  transformation,  f*(^,  T))  •  This  method  has  been  modified  in  order  to 
improve  computational  efficiency.  The  results  are  presented  below  in  two 
parts:  the  first  part  treates  the  confutation  of  the  coordinates  (x  ,  y)  of  the 
subsatellite  point  (zero  scan  angle);  the  second  part  develops  the  necessary 
added  formalism  for  datum  points  lying  along  a  scan  arc,  as  defined  in  Sec¬ 
tion  3. 1.2.  Details  and  derivations  are  included  in  Appendix  C. 
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3.3. 1  Data  Base  Coordinates  of  the  Subsatellite  Point 


Computations  for  obtaining  the  data  base  coordinates  (x  ,  y)  of  the  subsatellite 
point  are  performed  as  follows: 

1.  livaluale  Equation  (3-1)  for  1  ,  m  ,  and  n  ,  and  note  the  appro¬ 
priate  cube  face  by  selecting  the  largest  of  1  ,  m  ,  and  n  and 
obsciwing  its  sign.  Tabic  3-2  establishes  the  necessary  conventions. 
The  general  variables  4*  »  t?*  »  Jtnd  p  are  the  geocentric  rectan¬ 
gular  coordinates  of  the  subsatellite  point  deduced  from  1  ,  m  , 

and  n  by  entering  Table  3-2,  where  p  is  the  largest  in  absolute 
magnitude  of  1  ,  m  ,  and  n  . 

2.  Compute  the  variables  a  and  j3  defined  by: 


a  =  4  Vp 


P  =  T?'/p 


(3-40) 


Table  3-2.  Cube  Face  Identification 


Face 

Number 

1 

77' 

P 

1 

m 

n 

1 

2 

-m 

n 

-1 

3 

-1 

n 

m 

4 

1 

n 

-m 

5 

-1 

-m 

n 

6 

-1 

m 

-n 
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3. 


Compute  X  and  y  from 


X  =  f*(a,  p) 
y  =  a) 


(3-41) 


where  £*(a  ,  /3)  is  developed  in  Appendix  C.  For  convenience,  the 
functional  relationship  is  repeated  below. 


f*(a,  /3)  =  r  <ay*  +  (1  -  y*)  +a^  (l~  a^) 


r  +  (r-M)(l-a^  +  (l-/5^)  C 


UO 

Ja^O 


(3-42) 


(1  -  a^)  0  +  (1  -  a^)  d.  a^*l 

uo  ‘  J 


where  Ot  and  P  are  defined  1^  Equations  (3-40),  =  1>/^,  and 

the  remaining  constants  are  defined  in  Appendix  C  (Equations  (C-10) 
and  (C-11)). 

3.3.2  Data  Base  Coordinates  for  Nonzero  Scan  Angle 

The  data  base  coordinates  (x  ,  y)  for  any  pair  of  datum  pointa  lying  on  a  scan 
arc  at  an  angular  separation  of  iti>/2  from  the  subsatellite  point  are  obtained 
by  performing  the  following  computational  steps; 

1.  Obtain  the  geocentric  longitude,  8  ,  and  colatitude,  0  ,  of  the  sub¬ 
satellite  point  using  Equations  (3-1)  and  (3-3).  (If  the  ground  trace 
has  already  been  determined,  its  direction  cosines  will  already 
have  been  computed  from  Equations  (3-1). ) 
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2.  Compute  the  sines  and  cosines  of  the  geocentric  longitudes,  0’  , 
and  colatitudcs,  0'  ,  of  the  datum  points,  using  Equations  (3-4), 
(3-5),  and  (3-6),  and  the  results  of  step  1. 

3.  Using  the  sines  and  cosines  of  6*  and  0’  in  Equations  (3-3),  ob¬ 
tain  the  two  sets  of  direction  cosines  (1*  ,  m*  ,  and  n')  of  the 
datum  points. 

4.  Identify  the  appropriate  cube  face(s)  and  variables  4'  •  P 

using  Table  3-2,  and  compute  values  for  a  and  using  Equa¬ 
tions  (3-40)  for  each  point. 

5.  Calculate  the  (x,  y)  coordinates  using  Equations  (3-41). 

3.4  APPROXIMATE  FORMS 

The  transformations  outlined  in  the  previous  subsections  share  the  common 
•  disadvantage  of  being  complex  and  difficult  to  evaltiate.  While  studying  the 
transformations  f  and  f*  for  the  previous  report,  it  was  determined  that  a 
linear  interpolation  scheme  was  sufficient  to  find  points  between  known  bench¬ 
mark  grid  points  (Reference  1,  Section  4.3).  It  is  therefore  evident  that  the 
transformation  equation  may  be  approximated  by  low-order  polynomials. 
subsection  presents  a  method  for  accurately  approximating  the  transformations 
over  the  entire  cube  face  using  polynomials  of  varying  degree  with  locally  eval¬ 
uated  and  tabulated  coefficients. 

Because  of  symmetry  properties  of  the  transformations,  it  is  necessary  to  eval¬ 
uate  the  coefficients  of  the  approximation  polynomial  over  only  one  octant  of 
one  face.  In  addition,  since 


x  =  f  y  =  f  (T?.^) 
C  =  f*(Xty);  T)=  f*(y,x) 


(3-43) 
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only  one  set  of  coefficients  is  necessary  to  determine  both  x  and  y  ,  and  a 
second  set  to  determine  ^  rj . 

From  Table  4-2  of  Reference  1,  it  can  be  seen  that  a  double  linear  interpola¬ 
tion  scheme  will  allow  computation  of  intermediate  points  to  an  accuracy  of 
0.25  pixel  over  a  range  of  at  least  16  pixels  (and  up  to  1024  pixels  in  the  central 
region)  on  a  4096  x  4096  square.  Consequently,  it  is  reasonable  to  assume  that 
a  higher  order  polynomial  would  allow  one-pixel  accuracy  over  considerably 
greater  distances.  As  a  first  approximation,  it  is  assumed  that  determination 
and  tabulation  of  coefficients  on  a  grid  with  64-pixel  intervals  will  suffice.  This 
interval  is  convenient,  since  it  corresponds  to  the  proposed  standard  data  base 
record  size. 

To  obtain  suitable  interpolating  polynomials,  the  standard  bivariate  interpola¬ 
tion  formulae  were  examined  for  three,  four,  and  six  points  (Reference  4, 
Sections  25. 2. 65  to  25.2. 67).  From  Figure  3-10,  the  formulae  are 

•  Three  Points  (Linear) 

'a  ^0  *  ■  P  ■  '>*  'o.O  *  '■'l.O  *  “'o.  I  * 

•  Four  Points 

‘a  <*0  *  =  0  -  P>  <1  -  1)  *  P(l  -  q) 

i 

+  q(l  -  p)  f^^  ^  +  pqf^^  j  +  O(h^)  (3-45) 
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Figure  3-10.  Grid  Points  Used  for  Interpolation 


I* 
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six  Points 


f  (Xn  Ph,  y  +  qk)  =  f  +  ~  f 

a  '  0  ’  *"0  ^  '  2  0,-1  2  -1,0 


+  (1  +  pq  -  -q^ 

^  p(p  -  2q  +  1)  . 

2  1,0 


(3-46) 


+  q(q  -  2p  +  1)  j  +  pqf^^  ^  +  0(h®) 


Equations  (3-44)  and  (3-45)  can  be  written  as  follows: 


fa  '  “w  *  *  “Ol  *  “OO 

(3-47) 

t,  (x.y)  =  'V  *  “lo  *  -  »„1  y  *  >00 

(3-48) 

The  coefficients  a^^  ,  a^^  ,  . . . ,  are  evaluated  by  setting  up  a  set  of  linear 
equations  corresponding  to  the  interpolation  benchmark  points  and  solving  them 
simultaneously. 

The  six-point  formula  is  less  convenient  to  apply  systematically,  since  it  re¬ 
quires  knowledge  of  points  outside  its  range  for  solution.  A  compromise  form 


'a 


20 


X  + 


y 

02 


(3-49) 


can  be  evaluated  on  five  points,  as  shown  in  Figure  3-11.  This  form  eliminates 
the  xy  cross-product  term  and  approximates  the  functional  form  with  a  parabo¬ 
loid  of  revolution  whose  axis  is  perpendicular  to  the  xy  plane. 
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1 


Figure  3-11. 


Five-Point  Solution 


A  computer  program  was  written  to  solve  the  system  of  simultaneous  equations 
arising  from  each  of  the  above  approximate  fonns  on  64  x  64  pixel  intervals  on 
one  octant  of  the  cube  face  (Figure  3-12).  A  standard  IBM  Scientific  Subroutine 
Package  (Reference  5)  routine,  SIMQ^,  was  used  for  solution  of  the  equations. 
To  test  the  accuracy  of  each  approximation,  a  set  of  evaluation  points  was 
chosen  (Figure  3-13).  The  approximations  were  used  successively  to  deter¬ 
mine  the  coordinates  of  each  evaluation  point  and  a  residual  was  computed  for 
each  point  by  solving  the  full  transformation.  Any  64  x  64  region  in  which  the 
largest  residual  exceeded  one  pixel  was  flagged  and  the  next  higher  degree  ap¬ 
proximation  was  tried.  Table  3-3  shows  the  number  of  coefficients  necessary 
for  one  pixel  accuracy  as  a  function  of  location. 

The  table  size  required  is  not  excessive  because,  as  stated  above,  only  one 
octant  (plus  diagonal  elements)  must  be  stored.  For  a  64  x  64  grid,  this  would 
be 


M  X  (31^ 


+  16)  =  M  X  977 


1 


SIMQ  employs  the  Gaussian  elimination  method  using  largest  pivotal  divisors. 
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Figure  3-12.  Origins  of  Coefficient  Set  Domains 
Necessary  To  Define  Entire  Face 


Table  3-3,  Number  of  Coefficients  Nece8sar3^  for  One-Pixel 
Accuracy  as  a  Function  of  Position 


3 

3  3 
3  3  3 
3  3  3  3 
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33333333 
333333333 
3333333333 
33333333333 

2048,  2048 
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4 
4  4 
3  4  4 
3  3  4  4 
3  3  3  3  4 
3  3  3  3  4 
3  3  3  3  3 
3  3  3  3  3 
3  3  3  3  3 
3  3  3  3  3 
3  3  3  3  3 
3  3  3  3  3 
3  3  3  3  3 
3  3  3  3  3 
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4096,  2048 
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where  M  is  the  largest  number  of  coefficients  required.  For  the  case  where 
M  =  5  ,  this  requires  4885  words. 

The  program  which  evaluates  the  coefflcients  for  the  approximate  direct  trans¬ 
formation  f  ^  is  showTi  in  Appendix  E.  Appendix  D  tabulates  the  coefficients 
derived  and  indicates  the  necessary  number  of  coefficients  for  one -pixel  ac¬ 
curacy.  The  use  of  these  approximations  will  be  discussed  in  Section  4. 1. 


X 


BENCHMARK 

POINTS 

EVALUATION 

POINTS 


X 


X 


Figure  3-13.  Arrangement  of  Evaluation  Points 


3.  5  COMPARISON  OF  TRANSFORMATIONS 

In  the  previous  subsections,  four  different  formulations  were  developed  for  the 
f*  function  (Inverse  transformation): 

1.  The  original  Cartesian  form  (Equation  (2-5)) 

2.  The  revised  Cartesian  form  (Equation  (3-42)) 

3.  The  polar  form  (Section  3. 2) 

4.  The  approximate  form  (Section  3. 4) 
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In  choosing  between  forms  1  and  2,  v.Jjich  produce  the  same  mapping,  form  2 
is  clearly  arithmetically  simpler  and  has  proved  to  produce  a  more  accurate 
determination  of  its  coefficients. 

Form  4  may  be  used  to  approximate  the  Cartesian  or  the  polar  forms  and  is  the 
most  efficient  form  computationally.  It  does,  however,  require  core  storage 
for  the  coefficient  tables. 

Figures  3-14  and  3-15  illustrate  the  differences  in  the  nature  of  the  mappings 
produced  by  the  Cartesian  and  polar  forms.  The  salient  difference  is  that  the 
polar  form  exhibits  nondifferentiability  at  crossings  of  the  diagonals  of  the 
face.  This  directional  discontinuity  arises  from  the  condition  that  all  great 
circles  pass  through  the  center  of  a  face  map  as  straight  lines.  This  problem, 
in  conjunction  with  the  somewhat  slower  computation  time  shown  in  Table  3-4, 
favors  the  Cartesian  form  for  most  applications.  The  polar  form  retains  the 
advantage  of  exactness  in  that  it  produces  equal  area  to  any  desired  accuracy 


Figure  3-14.  Distortions  Arising  From  the  Cartesian  Form 
of  the  Transformation 
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Figure  3-15.  Distortions  Arising  From  the  Polar  Form  of 
the  Transformation 


Table  3-4.  Execution  Times  for  Transformation 
Forms  on  the  IBM  S/360-91  (One 
Ordered  Pair) 


Execution  Time 

Transformation 

(microseconds) 

Cartesian 

300 

Polar 

1000 

Approximate 

40 

without  recomputing  the  coefficients.  It  is  difficult  to  compute  coefficients  that 
produce  Cartesian  transformations  which  will  match  retransformed  points  below 
the  2-nautical-mile  level. 

Both  the  Cartesian  and  polar  formulations  have  their  own  uses:  the  Cartesian 
system  at  moderate  resolutions  and  the  polar  system  at  the  highest  resolutions 
over  small  areas  not  including  a  diagonal.  For  benchmark  transformations, 
the  speed  of  either  transformation  is  adequate,  given  a  reasonably  great  inter¬ 
polation  distance.  The  high  computational  efficiency  of  the  approximate  forms 
tends  to  favor  this  approach  for  direct  computation  of  data  base  coordinates  (or 
for  computation  of  corresponding  points  in  the  scan  line  buffer).  This  technique 
could  also  be  particularly  useful  in  the  comer  regions  where  Interpolation  is 
not  accurate  over  useful  distances. 
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SECTION  4  -  I/O  SYSTEM  DESIGN  CONSIDERATIONS 


4.1  ALTERNATIVE  METHOD  OF  FAST- FILLING 

Section  2.3  briefly  outlines  a  proposed  alternath’e  fast-filling  scheme  for  the 
preprocessor  program.  This  method,  described  in  Reference  3,  partially 
solves  the  problem  of  edge  and  vertex  crossings  and  eliminates  the  necessity 
for  explicit  data  spreading.  By  combining  this  approach  with  the  approximation 
techniques  descussed  in  this  report,  the  overall  efficiency  of  the  preprocessor 
function  can  be  increased  still  further. 

Sections  3. 1  and  3.3  present  a  method  for  obtaining  data  base  coordinates 
(X  ,  y)  directly  from  Keplerlan  orbital  elements  and  scanner  angle.  By  omit¬ 
ting  the  last  step  of  the  process  outlined  in  Section  3.3,  positions  in  the  Ifx  ,  j3) 
system  can  be  obtained,  which  is  simply  a  rescaled  (4  ,  t))  system  (ot  =  , 

P  Tj/r^)  .  Using  this  algorithm,  a  benchmark  grid  in  (oc  ,  j3)  can  be  obtained 
and  a  procedure  can  then  be  followed  which  consists  of  a  modification  to  that 
described  in  Reference  3.  As  outlined  in  Section  2.3  of  Reference  3,  filling  is 
performed  only  for  complete  data  base  records  (thus  avoiding  the  necessity  to 

read  and  update  partially  covered  records).  The  procedure  employs  a  new 

R  R  R 

rotating  Earth-centered  coordinate  system  (x  ,  y  ,  z  )  in  which  both  scan¬ 
ner  and  data  base  benchmarks  are  defined.  CSC  proposes  modifying  this  pro¬ 
cedure  as  follows: 

R  R  R 

1.  Instead  of  using  the  x  ,  y  ,  z  system,  both  scanner  and  bench¬ 
mark  grids  in  a  coplanar  system  should  be  defined,  projecting  the 
scan  line  grid  into  the  inscribed  cube  face  and  transforming  the  data 
base  coordinates  to  the  same  (ct  ,  P)  system. 

2.  The  data  base  benchmarks  should  be  defined  in  terms  of  the  pro¬ 
jected  scanner  benchmarks  (i.e. ,  the  (E,  L)  coordinates  of  the 
data  base  benchmarks  should  be  defined  as  shown  in  Figure  4-1). 


4-1 


n 


•  DATA  BASE  BENCHMARKS 

X  SCANNER  BENCHMARKS 

a.  fi  CUBE  FACE  COORDINATES 

E,  L  SCANNER  COORDINATES 

X,  y  DATA  BASE  COORDINATES 


Figure  4-1.  Coplanar  Representation  of  Scanner  and  Data  Base  Grids 


3.  The  (E,  L)  coordinates  of  each  data  base  location  should  be  inter¬ 
polated  and  the  nearest  pixel  stored  in  that  location.  (This  may  be 
done  by  offsetting  the  (E,  L)  system  1/2  pixel  in  both  coordinates 
and  obtaining  E  and  L  indices  of  the  desired  pixel  by  truncation. ) 

The  objective  of  this  modification  is  to  avoid  the  necessity  of  searching  for  the 
nearest  pixel  in  a  three-dimensional  space. 

Neither  this  method  nor  the  procedure  in  Reference  3  completely  eliminates  the 
edge  and  vertex  crossing  problems.  It  is  necessary  to  predetermine,  from  the 
benchmarks  of  the  scanner  array,  which  data  base  records  will  be  completely 
covered  by  the  scanner  data  currently  in  the  input  buffer.  This  process  must 
check  for  edge  and  vertex  crossings. 

Figure  4-2  illustrates  the  relationship  between  the  scanner  data  input  buffer 
and  the  data  base  records.  The  scan  buffer  must  have  sufficient  length  (in  the 
along-track  direction)  to  cover  the  longest  diagonal  of  a  data  base  record.  For 
the  64  X  64  pixel  case,  this  will  be  approximately  128  pixels  in  the  comers  of 
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the  face.  The  scan  line  buffer  must  be  organized  as  a  creeping  buffer,  and  no 
more  lines  may  be  added  in  the  forward  direction  (and  deleted  in  the  backward 
direction)  than  the  number  of  lines  by  which  the  buffer  exceeds  the  minimum 
needed  to  cover  the  longest  diagonal.  A  typical  scanner  input  buffer  would  have 
a  width  (along  the  scan  direction)  of  approximately  1600  nautical  miles  (about 
1344  pixels),  with  a  minimum  length  (in  the  ground  track  direction)  of  128  pixels. 
This  \yill  require  sufficient  buffer  storage  for  172,032  pixels.  Two  methods 
emerge  for  handling  this  large  array.  The  first  is  to  pack  more  than  one  pixel 
per  60-bit  word.  This  would  allow  the  entire  buffer  to  be  maintained  in  central 
memory  on  the  CYBER/175  but  would  introduce  a  large  packing  and  unpacking 
cost.  A  second  method  would  be  to  maintain  the  buffer  in  ECS,  segmented  so 
that  it  could  be  brought  into  core  In  blocks  as  required.  The  scheduling  of  these 
blocks,  along  with  the  ordering  of  output  records,  could  be  determined  from 


Figure  4-2.  Relationship  of  Scanner  Input  Buffer  and  Data  Base  Records 
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the  (X  ,  y)  locations  of  the  projected  and  transformed  scanner  input  benchmark 
grid,  ordering  the  filling  of  data  base  records  by  their  corresponding  along-scan 
locations. 

Because  the  objective  of  this  new  method  is  to  reduce  redundant  stores,  it 
could  also  be  useful  to  consider  windowing  down  the  along-scan  length  of  the  in¬ 
put  buffer  as  a  function  of  latitude  to  reduce  scan  path  overlap  at  high  latitudes. 
Similarly,  at  low  latitudes,  it  may  be  necessary  to  subdivide  data  base  records 
into  two  or  more  longitude  bands  and  partially  update  the  records  rather  than 
fully  rewriting  them. 

4.2  MAINTENANCE  OF  MULTIPLE-  AND  VARIABLE-RESOLUTION  DATA 
BASES 

Considerable  savings  in  computation  and  I/O  time  and  in  data  base  storage  can 
be  realized  by  reducing  the  resolution  of  the  data  base  in  areas  of  lesser  inter¬ 
est  (such  as,  perhaps,  over  land  and  in  the  southern  hemisphere).  It  may  also 
be  desirable'to  be  able  to  maintain  higher  resolution  lor  other  selected  areas. 
One  method  for  providing  such  variable  resolution  coverage  consists  of  main¬ 
taining  a  single,  fully  global  data  base  at  the  lowest  desired  resolution,  plus 
separate  higher  resolution  partial  data  bases.  A  second  method  consists  of 
maintaining  a  single  data  base  of  variable  resolution. 

In  either  case,  it  will  be  necessary  to  provide  a  system  for  indicating  the  de¬ 
sired  resolution  as  a  function  of  global  position.  This  can  be  provided  for  by 
establishing  a  data  block  properties  table  (for  example  at  the  128  x  128  or 
256  X  256  element  level)  which  would  contain  the  desired  resolution  described 
as  a  multiple  Of  the  finest  usable  resolution  or  by  the  number  of  bits  required 
for  the  serial  address.  This  would  allow  the  preprocessor  program  to  deter¬ 
mine  the  desired  resolution  at  any  point. 

For  the  multiple  data  base  scheme,  the  preprocessor  program  would  first  de¬ 
termine  which  elements  of  the  resolution  table  are  covered  by  the  data  currently 
in  the  scan  line  buffer  and  then  access  the  table  to  determine  what  actions  must 
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be  taken.  This  table  would  consist  of  entries  for  each  of  the  desired  resolution 
levels  indicating 

•  The  desired  resolution 

•  A  base  address  for  the  block  of  records  corresponding  to  the  table 
element 

•  The  size  of  the  data  base  records  (in  elements) 

The  resolution  and  record  size  for  ail  records  in  a  block  corresponding  to  one 
table  element  can  both  be  expressed  as  a  number  of  bits  to  be  added  to  the 
corresponding  binary  division  levels  of  the  basic  global  data  base,  i.  e. ,  if  the 
desired  resolution  is  twice  the  linear  resolution  of  the  basic  data  base,  two 
more  address  bits  must  be  added  and,  if  the  records  are  to  have  the  same  num¬ 
ber  of  elements,  the  record/element  division  of  the  serial  string  must  take 
place  two  bits  lower. 

The  serial  string  for  each  resolution  level  would  be  kept  separate,  and  the  use 
of  a  base  address  for  each  record  block  would  obviate  any  need  for  unused  filler 
records  in  the  data  base.  For  maximum  flexibility,  each  resolution  level  could 
be  maintained  as  a  separate  file. 

With  multiple  data  bases,  the  size  of  a  data  base  record  would  be  adjustable  to 
allow  an  optimum  match  of  record  size  to  scanner  iigiut  buffer  size  and  data 
item  size  (where  one  data  item  might  be  all  of  the  simultaneously  available, 
colocated  scanner  readings).  To  maintain  this  match,  a  separate  pass  through 
the  data  will  be  required  for  each  resolution  level. 

In  the  variable  resolution  approach,  only  one  pass  through  the  data  is  required 
and  all  data  records  cover  the  same  area.  This  means  that  the  ratio  of  the 
scanner  buffer  size  to  the  data  record  area  remains  fixed,  as  do  the  total  num¬ 
ber  of  I/O  operations.  The  total  number  of  data  records  remains  constant; 
thus,  a  consistent  record  numbering  and  ordering  can  be  maintained,  avoiding 
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the  necessity  of  a  table  lookup  for  a  block  address,  as  in  the  multiple  data  base 
scheme.  In  this  case,  the  data  block  lookup  table  need  contain  only  the  reso¬ 
lution  (which  will  also  imply  the  record  length). 

Both  of  the  above  schemes  have  advantages  and  disadvantages: 

•  The  multiple  data  base  scheme  requires  redundant  stores  and  com¬ 
putation 

•  The  variable  resolution  scheme  requires  stringing  together  unlike 
data  base  records,  thus  requiring  a  variable  record  length  capability 

•  The  multiple  data  base  scheme  allows  easy  addition  or  deletion  of 
nonstandard  resolution  areas  for  coverage  of  new  areas  of  interest 

•  The  variable  resolution  scheme  allows  more  systematic  control 
over  the  data  base  storage  allocations 

A  possible  resolution  of  these  conflicts  is  to  use  a  hybrid  system  in  which  the 
variable  resolution  format  would  be  used  for  resolutions  coarser  than  a  set 
standard  resolution  level,  and  the  multiple  data  base  approach  for  higher  reso¬ 
lution  levels  (and  for  other  subdata  bases  required  for  special  purposes  or 
limited  durations). 

A  unified  block  properties  table  would  be  maintained  which  would  allow  display 
of  the  current  status  and  availability  of  a  data  record  in  the  system.  This  table 
could  be  expanded  to  include  such  information  as  the  last  update  time,  satellite 
ID,  dayside/nightslde  information,  etc.  The  advantages  of  such  a  scheme  would 
be  that  multiple  passes  through  the  data  would  be  restricted  to  selected  regions 
of  highest  resolution,  special  changes  in  data  base  content  could  be  accomplished 
without  restructuring  the  main  data  base  as  a  whole,  systematic  control  over 
the  entire  data  base  could  be  maintained,  and  the  serial  string  addressing 
scheme  would  be  maintained  for  the  standard  resolution  data. 
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SECTION  5  -  GEOGRAPHIC  ORIENTATION 


5. 1  INTRODUCTION 

In  order  to  maximize  the  use&lness  and  efficiency  of  the  spherical  cube  data 
base,  it  may  be  necessary  to  orient  the  inscribed  cube  in  some  arbitrary  man¬ 
ner  witli  respect  to  the  Earth.  The  simplest  orientation  of  the  cube  aligns  one 
Cartesian  axis  of  the  cube  coordinate  system  with  the  Earth's  axis  of  rotation 
and  points  another  cube  axis  toward  the  prime  meridian.  A  simple  variation 
is  to  have  one  cube  axis  aligned  with  the  Earth's  axis  of  rotation  and  another 
pointing  to  a  specified  longitude.  This  section  will  discuss  these  orientations 
and  will  present  an  algorithm  for  producing  arbitrary  cube  rotations. 

5. 2  SIMPLE  ORIENTATIONS 

In  Section  4. 2  and  Appendix  E.  6  of  Reference  1,  algorithms  are  presented  for 
the  transformation  of  Earth-centered  Cartesian  coordinates  (a  ,  b  ,  c)  to  data 
base  coordinates  (x  ,  y)  .  For  the  case  where  the  Earth-centered  (a  ,  b  ,  c) 
system  is  aligned  as  shown  in  Figure  5-1,  the  Cartesian  coordinates  can  be 
represented  as: 


a  =  cos  (latitude)  cos  (longitude) 

(6-1) 

b  =  cos  (latitude)  sin  (longitude) 

(5-2) 

c  =  sin  (latitude) 

(5-3) 

To  accomplish  a  simple  rotation  about  the  Earth's  polar  axis,  the  following 
algorithms  can  be  used: 

a  =  cos  (latitude)  cos  (longitude  +  6)  (6-4) 
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b  cos  (latitude)  sin  (longitude  +  6) 


(5-5) 


c  =  sin  (latitude)  (5-6) 

where  6  represents  the  angle  of  rotation. 


N  (POLAR  AXIS) 


LONGITUDE  SO**E 


Figure  5-1.  Orientation  of  the  Axes  of  the  Cube  (Unrotated) 


5.3  THE  GENERAL  ROTATION  CASE 

In  order  to  conveniently  generate  arbitrary  rotations  of  the  data  base  cube  with 
respect  to  the  sphere,  a  method  was  adopted  for  making  successive  rotations 
of  the  sphere  about  each  of  the  axes  of  the  cube.  By  choosing  the  order  in  a 
suitable  manner,  this  yields  a  conceptually  simple  way  of  arriving  at  a  given 
orientation. 
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Taking  a  general  set  of  orthogonal  axes,  G  ,  G  ,  G  ,  as  shown  in  Figure  5-2, 

J,  A  U 

the  rotation  about  axis  G^  can  be  expressed  as  follows: 

Let  a  be  the  desired  angle  of  rotation  about  G^^  .  For  any  arbitrary  point, 
g,  •  gn  •  go  >  position  may  be  expressed  in  terms  of  the  following  three 

X  '  A  ‘  i5 

variables : 

•  gj^  ,  the  Gj^  coordinate 

•  ^  ,  the  angular  position  of  the  projection  of  g^^  ,  g^  t  g^  on  the 
G  G  plane  measure  from  G  in  the  clockwise  sense 

A  u  O 

•  p  ,  the  perpendicular  distance  from  g^  ,  g2  •  g^  to  the  G^  axis 

For  a  pure  rotation  about  G^  (Figure  5-3),  g^  and  p  are  constant  and  the 

final  orientation,  g'  ,  gl  ,  gl  ,  is  determined  the  new  orientation  angle,  y  , 

X  A  U 

where 


y  =  a  +  j8 


and  since 


then 


y  =  a  +  tan"^  <83/82^ 


where 


P  =  g2+83 


(5-7) 


(5-8) 


(5-9) 


(5-10) 
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Figure  5-2.  Point  g]^  .  g2  i  gs  to  the  General  Reference 
Frame,  Gp  ,  G2  »  G3  . 


Figure  5-3.  Rotation  of  the  Point  (g^  ,  g2  •  ga) 
About  the  G^  Axis 
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yielding 


gi=gl 

(5-11) 

g^=pcos  y 

(5-12) 

g^=p8iny 

(5-13) 

Using  this  general  approach,  a  rotation  through  an  angle  is  first  made  about 
the  axis  coincident  with  the  Earth's  axis  of  rotation.  This  is  functionally  equiv¬ 
alent  to  the  rotation  described  in  Section  5. 1  and  results  in  a  rotation  of  the 
longitude  line  into  the  center  of  face  number  1. 

The  second  rotation  is  made  about  the  axis  normal  to  the  pole  and  emerging 
throus^  a  face  normal  to  face  number  1  (axis  b  in  Figure  5-1).  If  this  rotation 
is  through  an  angle  a,,  •  then  the  point  with  longitude  a  and  latitude  a 
will  be  rotated  into  the  center  of  face  1. 

A  third  rotation  about  the  axis  emerging  from  the  center  of  face  1  will  produce 
a  "position  angle"  change  on  face  1  but  will  not  translate  the  point  at  the  center 
of  the  face. 

Using  this  method,  any  point  whose  latitude  and  longitude  are  known  may  be 
translated  to  the  center  of  face  1  and  the  rotational  aspect  of  that  face  may  be 
altered  without  changing  the  latitude  and  longitude  of  the  center.  Referring  to 
Figure  5-1,  the  order  of  these  rotations  is  about  axes  c,  b,  and  a,  respectively. 

Appendix  E  contains  an  algorithm  for  the  application  of  the  general  transforma¬ 
tion  to  successive  axes  through  the  interchange  of  components. 

Figures  5-4  and  5-5  illustrate  the  continental  outlines  after  the  indicated  rota¬ 
tions  and  are  typical  of  maps  generated  by  this  process. 
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Figure  5-4.  Data  Base  Map  of  Earth  With  Cube  Orioited  90  Degrees 
West  of  Greenwich 
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5-5.  Data  Base  Map  of  Earth  With  Cube  Oriented  90  Degrees 
West  of  Greenwich  and  30  Degrees  South  of  the  Equator 


SECTION  6  -  SUMMARY  AND  CONCLUSIONS 


6.1  NEW  RESULTS 

The  purpose  of  the  research  reported  herein  is  to  improve  the  procedures  and 
algorithms  required  for  the  efficient  maintenance  of  data  in  a  quadrilateralized 
spherical  cube  Earth  data  base.  The  new  results  obtained  relate  to  the  fast¬ 
filling  of  the  data  base.  A  technique  was  developed  for  rapidly  finding  the  groimd 
trace  and  scanner-sensed  locations  from  Keplerian  elements  and  scanner  angle. 
The  previously  derived  form  of  the  inverse  transformation,  f*  ,  was  improved 
by  a  change  of  variables.  A  new  inverse  transformation  was  derived  using  polar 
coordinates,  and  a  technique  was  developed  for  approximating  the  transforma¬ 
tion  fvmctions  using  low-order  polynomials. 

All  of  these  results  will  contribute  to  improving  the  accuracy  and  speed  of  the 
preprocessor  function.  The  general  functioning  of  the  preprocessor  input/output 
algorithms  was  re-examined  in  light  of  an  alternative  scheme  (Reference  3). 

This  scheme  and  the  previously  developed  fast-filling  method  (Reference  1) 

were  compared,  and  a  modification  of  the  scheme  described  in  Reference  3  is 

! 

recommended.  ; 

i 

In  addition,  improved  world  maps  in  data  base  coordinates  were  produced  fo:b 
use  in  determining  an  optimum  orientation  of  the  data  base  cube.  Additional 
efforts  \s  ere  devoted  to  improving  the  accuracy  of  the  transformations  described 
in  Reference  1.  The  revised  coefficients  for  these  transformations  are  included 
in  Appendix  B  of  this  report. 

6.2  RECOMMENDATIONS  FOR  FURTHER  STUDY 

A  number  of  topics  considered  during  this  study  are  worthy  of  additional  atten¬ 
tion.  The  polar  formulations  of  the  transformation  functions  should  be  further 
developed  using  other  constraints  to  maximize  conformality  of  the  mapping. 


6-1 


The  role  of  the  MAP  HI  processor  in  Implementing  the  preprocessor  function 
was  excluded  from  this  study  because  documentation  is  not  available.  The  role 
which  MAP  in  and  the  high-speed  communications  system's  processor  could 
assume  in  evaluating  the  transformation  equations,  particularly  the  approximate 
transformation,  shovild  be  explored  further. 

The  various  CYBER  I/O  access  methods  should  be  studied  further,  particularly 
to  ascertain  whether  or  not  they  may  cause  fragmentation  of  the  data  base. 

Consideration  should  be  given  to  producing  display  output  in  data  base  coordi¬ 
nates  with  a  geographic  grid  overlay,  thereby  saving  a  major  computation  cost 
in  retransforming  the  data. 
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APPENDIX  A  -  GROUND  TRACE  EQUATIONS 


This  appendix  presents  the  derivation  of  the  expressions  summarized  in  Sec¬ 
tion  3. 1  and  referred  to  elsewhere  in  this  document.  Two  problems  are  ad¬ 
dressed:  (1)  the  computation  of  the  geocentric  direction  cosines  (1 ,  m  ,  n)  of 
the  subsatellite  point  in  terms  of  the  Keplerian  orbital  elements  of  the  satellite 
and  (2)  the  computation  of  the  geocentric  longitudes  and  colatitudes  of  any  pair 
of  datum  points  on  a  scan  arc  in  terms  of  previously  conq>uted  quantities  (e.  g. , 

1  ,  m  ,  n)  .  A  spherical  Earth  is  assumed  throughout. 

A.  1  THE  SUBSATELLITE  POINT 

Figure  A-1  shows  the  geometric  relationship  of  the  satellite,  the  satellite's 
orbit,  and  the  Earth's  equator  to  the  (rotating)  rectangular  coordinate  system 
(X  ,  y  ,  z)  .  The  coordinate  system  is  fixed  in  the  Earth  and,  therefore,  ro¬ 
tates  at  the  diurnal  rate.  The  K^lerian  orbital  elements  whose  definitions 
are  given  by  the  figure  are  the  inclination,  i  ;  the  argument  of  perigee,  to  , 
shown  as  part  of  the  angle  (co  +  f)  ,  where  f  is  the  satellite's  true  anomaly; 
and  the  node,  Q  ,  which  is  not  the  usual  fixed  quantity,  but  varies  with  the 
Fldereal  time  (diurnal  rate),  s  ,  as  follows: 

0  =  0  -  s  (A-1) 

o  ' 

where  O^  is  the  value  of  0  at  an  arbitrary  initial  time. 

Figure  A-1  shows  the  subsatellite  point  having  geocentric  coordinates  longitude, 
8  ,  and  colatitude,  0  ,  and  lying,  by  definition,  on  the  line  joining  the  satellite 
and  the  Earth's  center.  Therefore,  the  direction  cosines  (1  ,  m  ,  and  n)  de¬ 
scribing  the  direction  to  the  satellite  in  the  x,  y,  z  system  are  identical  with 
those  for  the  subsatellite  point.  The  straightforward  relationship  between  the 
direction  cosines  and  6  and  0  will  be  e}q;)lolted  in  Section  A.  2.  First,  the 
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Figrure  A-1.  Satellite  Orbit  Orientation  and 
Geographic  Coordinates 
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expressions  for  1  ,  m  ,  and  n  in  terms  of  O  ,  co  .  i  ,  and  f  are  derived 
below. 

Consider  the  geocentric  position  vector,  ■?,  of  the  satellite  in  its  orbit.  Ex¬ 
press  T  in  the  two-dimensional  rectangular  coordinate  system  (p  ,  q)  in  the 
orbital  plane,  where  the  unit  vector  p  is  in  the  direction  of  the  node  and  the 
unit  vector  q  is  in  the  direction  90  degrees  in  advance  of  p  (see  Figure  A-1) 

■f  =  r  cos  (oj  +  f)  p  +  r  sin  (0)  +  f)  q 

where  r  is  the  magnitude  of  "r  .  The  coordinates  (r^  ,  r^  ,  r^)  of  the  satel¬ 
lite  in  the  geocentric  equatorial  rectangular  coordinate  system  are.  then, 

r  =?  •  X  =  r  cos  (w  +  f)  p  •  x  +  r  sin  (o)  +  f)  q  •  x 

X 

r^  =■?•  y  =  r  cos  (o)  +  f)  p  •  y  +  r  sin  (co  +  f)  q  •  jr  (A-2) 

r^  =■?  ‘Z*  +rsin(UJ  +  f)q»z 

where  the  term  in  p  •  z  for  r  is  omitted,  since  p  •  “z  =  0  . 

z 

Now, 


p  •  X  =  cos  O 

q  •  X  =  cos  (90°  +  0)  cos  i  =  -sin  0  cos  i 
p  •  y  =  sin  O  (A-3) 

q  •  y  =  sin  (90°  +  0)  cos  i  =  cos  0  cos  i 
q  •  "z  =  cos  (90°  -  1)  =  sin  i 
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Substituting  Equations  (A-3)  into  Equations  (A--2),  ,  and  r 

r  =  r  cos  O  cos  (w  +  f)  -  r  sin  O  sin  (co  +  f)  cos  i 

X 

r  =  r  sin  0  cos  (o)  +  f)  +  r  cos  O  sin  (to  +  f)  cos  i 

y 

=  r  sin  (to  +  f)  sin  i 

Since  the  satellite  position  vector,  r" ,  in  the  y,  z)  system  is 


r  =  r  x  +  ry  +  r  z 
X  jr  z 

the  direction  cosines  (1 ,  m  ,  n)  of  its  position  are  defined  by 


r‘.x=  1^1  lx  11 
^•y  =  1^1  |y  I  m 


Now,  from  Equation  (A-5), 


■  X 


=  r  X 
X 


X 


r 

X 


r-y  =  y  .y  =  ry 

/S  /% 

r»z  =  r  z*z  =  r 
z  z 


A-4 


become 


(A-4) 


(A-5) 


(A-6) 


which,  when  substituted  into  Equations  (A-6),  yields 


l  =  |xl  =  ryr 

m  =  ty/I-ri  ly  I  =  ry/r  (A-7) 

n  =  r  /|T|  |z  j  =  r  /r 

since  |x  1  =  |y  1  =  |z  |  =  1  ,  and  |r“|  =  r. 

Equations  (A-4)  and  (A-7)  then  yield 

1  =  cos  O  cos  (CO  +  f)  -  sin  0  sin  (u)  +  f)  cos  i 
m  =  sin  O  cos  (co  +  f)  +  cos  O  sin  (co  +  f)  cos  i  (A-8) 

n  =  sin  (CO  +  f)  sin  i 

Equations  (A-8)  provide  one  way  of  defining  the  "ground  trace"  of  the  satellite, 
since  1  ,  m  ,  and  n  mry  be  regarded  as  the  "unit"  geocentric  coordinates  of 
the  subsatellite  point.  That  is,  assuming  a  spherical  Earth,  the  actual  coordi¬ 
nates  of  the  subsatellite  point  are  obtained  by  multiplying  1  ,  m  ,  and  n  by 
the  Earth's  radius  in  whatever  units  are  required  (e.g. ,  kilometers).  Equa¬ 
tion  (A-8)  is  given  as  Equation  (3-1)  in  the  text. 

A.  2  THE  LOCATION  OF  DATUM  POINTS  ON  A  SCAN  ARC 

The  satellite's  onboard  sensors  are  assumed  to  sweep  out  a  scan  arc  having 
ground  trace  geometric  properties  as  defined  in  Section  3. 1. 2  and  as  illustrated 
in  Figure  A-2.  It  is  convenient  to  regard  the  ground  trace  of  the  scan  arc  as 
being  composed  of  pairs  of  datum  points  symmetrically  placed  on  either  side 
of  the  subsatellite  point,  and  to  regard  (|)/2  as  a  general  variable  locating  the 
points  by  assigning  alternate  algebraic  signs  to  (/)/2  (e.  g. ,  the  ends  of  the 
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QROUNO  TRACE 


Figure  A-2.  The  Scan  Arc  as  a  Central  Angle 
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scan  arc  are  located  at  from  the  subsatellite  point).  It  is  desired  to 

obtain  pairs  of  values  for  the  geocentric  longitude,  6'  ,  and  colatitude,  0'  ,  of 
the  datum  points  in  terms  of  known  quantities  in  a  computationally  efficient 
manner. 

The  sines  and  cosines  of  the  geocentric  longitude,  6  ,  and  colatitude,  0  ,  of  the 
subsatellite  point  are  required  as  intermediate  quantities  in  terms  of  1  ,  m  , 
and  n  from  Equations  (A- 8).  By  inspection  of  Figure  A-1, 

cos  0  sin  0  =  1 

sin  6  sin  0  =  m  (A-9) 

cos  0  =  n 

Care  must  be  taken  in  using  Equations  (A>9)  to  obtain  sin  0  ,  cos  0  ,  sin  6  , 
and  cos  0  ,  since  the  algebraic  signs  of  these  quantities  must  be  preserved  for 
later  use.  This  is  readily  achieved  by  noting  that  cos  0  yields  0  uniquely, 
since  0^0^  180  degrees  .  Straightforward  manipulation  yields  the  remaining 
quantities.  Equation  (A-9)  is  given  in  the  text  as  Equation  (3-3). 

To  obtain  the  values  of  0’  ,  the  necessary  relations  may  be  derived  by  con¬ 
sidering  the  two  spherical  triangles  PSQ  and  PQR  shown  in  Figure  A-2  and 
having  the  common  side  PQ  ,  which  is  the  colatitude,  0  ,  of  the  subsatellite 
point.  These  triangles  are  shown  in  detail  in  Figure  A-3,  where  it  is  noted 
that  the  arc  PS  is  90  degrees,  that  the  scan  arc  intersects  the  ground  trace 
at  right  angles,  and  that  point  R'  may  be  regarded  as  any  point  on  the  scan 
arc  at  a  distance  of  (|)/2  from  the  subsatellite  point,  Q  .  (Figure  A-2  shows 
only  the  end  point  of  the  scan  arc,  R  ,  at  a  distance  0^/2  from  the  subsatellite 
point. )  Also  shown  in  Figure  A-3  are  two  auxiliary  angles,  A  and  /3  ,  such 
that  A  +  /3  +  90  degrees  =  360  degrees  . 
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Figure  A-3.  Two  Adjoining  Spherical  Triangles  and 
the  Subsatellite  Point  Q 
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The  cosine  law  applied  to  triangle  PQR*  yields 


cos  0'  =  cos  0  cos  0/2  +  sin  0  sin  0/2  cos  /3  (A-10) 

Since  0/2  is  a  quantity  given  a  priori,  all  of  the  quantities  to  the  right  of  the 
equality  in  Equation  (A-10)  are  known  except  cos  /3  ,  which  is  obtained  as  fol¬ 
lows. 

By  definition 


/3  =  270°  -  A 


Therefore 

cos  =  cos  (270°  -  A)  =  cos  (270°)  cos  A  +  sin  (270°)  sin  A 

(A-11) 

=  -sin  A 


Applying  the  sine  law  to  triangle  PSQ  , 


sin  A 
sin  90° 


sin  (90°  -  i) 
sin  0 


or 


(A-12) 


sin  A  =  cos  i/sin  0 

Employing  Equations  (A-11)  and  (A-12)  in  Equation  (A-10)  yields 

cos  0'  =  cos  0  cos  0/2  -  cos  i  sin  0/2  (A-13) 
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which  is  the  desired  relation  for  cos  0*  .  Equation  (A-13)  yields  0'  viniquely, 
since  0  £  0'  £  180  degrees  ,  so  that  sin  0*  follows  immediately.  By  appropri¬ 
ately  redrawing  Figure  A-3  and  repeating  the  above  analysis,  it  can  be  readily 
seen  that  0*  for  the  datum  point  corresponding  to  -0/2  on  the  other  half  of  the 
scan  arc  is  obtained  from  Equation  (A-13)  simply  by  replacing  0/2  with  -0/2  . 

Next,  to  obtain  the  values  of  0'  ,  Figure  A-3  is  again  used  to  note  that  the 
cosine  law  applied  in  triangle  PQR'  yields 

cos  0/2  =  cos  0  cos  0'  +  sin  0  sin  0’  cos  (0'  -  0) 
or  (A-14) 

cos  (0*  -  0)  =  (cos  0/2  -  cos  0  cos  0')/8in  0  sin  0’ 

Since  0  and  all  other  quantities  in  Equation  (A-14)  are  previously  known 
throu^  Equations  (A-9),  it  remains  only  to  obtain  an  independent  relation  for 
sin  (0*  -  0)  to  uniquely  specify  the  quadrant  of  O'  .  The  sine  law  applied  to 
triangle  PQR'  yields 


sin  B  _  sin  (0'  -  6) 
sin  0*  sin  0/2 


or 


(A-15) 


sin  (0'  -  0)  =  sin  ^  sin  0/2/Bin  0' 


To  obtain  sin  ,  it  should  be  noted  that,  by  definition, 


Therefore, 


sin  /3  =  sin  (270°  -  A)  =  sin  (270°)  cos  A  -  cos  (270°)  sin  A 
=  -cos  A 


(A-16) 


Applying  the  cosine  law  for  angles  to  triangle  PQS  , 

cos  A  =  -cos  (90°  -  i)  cos  (0  -  0)  +  sin  (90°  -  i)  sin  (6  -  0)  cos  (90°) 

(A-17) 

=  -sin  i  cos  (6  -  O) 

However,  the  cosine  law  for  sides  in  triangle  PQS  yields 

cos  (w  +  f)  =  cos  (90°)  cos  0  +  sin  (90°)  sin  0  cos  (0  -  O) 
or 

cos  (0-0)  =  cos  (0)  +  f)/sln  0 

Therefore,  Equation  (A-17)  becomes 

cos  A  =  -sin  i  cos  (w  +  f)/sin  0  (A-18) 

Substituting  Equation  (A-18)  into  Equation  (A-16),  and  then  substituting  the 
result  into  Equation  (A-15)  yields 

sin  (0'  -  0)  =  cos  (w  +  f)  sin  i  sin  0/2/sin  0  sin  0'  (A-19) 


Equations  (A-14)  and  (A-19)  are  the  desired  expressions  that  uniquely  deter¬ 
mine  0'  from  previously  known  quantities.  Alteniative  forms  for  Equa¬ 
tion  (A-19)  were  found,  but  the  one  given  here  offers  the  simplest  functions 
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of  the  quantities,  all  of  which  are  previously  known  and  is  the  most  rapid  to 
compute,  since  no  additional  evaluations  of  trigonometric  functions  are  required. 

Again,  by  appropriately  redrawing  Figure  A-3  and  repeating  the  above  analysis, 
it  can  be  readily  seen  that  Equations  (A-14)  and  (A-19)  uniquely  yield  6'  for 
the  datum  point  corresponding  to  -^/2  on  the  other  half  of  the  scan  arc.  This 
is  accomplished  by  entering  Equations  (A-19)  with  -0/2  ,  and  entering  both 
Equations  (A-14)  and  (A-19)  with  the  appropriate  0*  from  Equation  (A-13). 

The  important  results  just  derived  are  Equations  (A-13),  (A-14)  and  (A-19), 
which  yield  0'  and  6'  for  pairs  of  datum  points  on  any  scan  arc,  and  are  cited 
as  Equations  (3-4),  (3-5),  (3-6),  respectively,  in  the  text.  Caution  is  neces¬ 
sary  when  applying  Equations  (A-14)  and  (A-19)  near  the  poles,  since  these 
relations  are  singular  for  subsatellite  points  or  datum  points  at  the  poles. 

This  is  expected,  since  the  geocentric  longitudes  6  and  6'  are  undefined  at 
the  poles.  No  problem  need  arise  in  computations,  however,  since  the  satel¬ 
lites  are  not  e:q)ected  to  have  exactly  90-degree  inclinations  and  since,  for  the 
case  of  an  equatorial  orientation  of  the  data  base  cube,  it  is  a  simple  matter 
to  regard  the  cube  face  coordinates  of  such  datum  points  as  zero  within  an 
arbitrary,  but  nonzero,  distance  from  the  poles. 
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APPENDIX  B  -  BEST  VALUES  FOR  THE  COEFFICIENTS  OF  f* 


For  the  form  of  the  inverse  mapping  function,  f*  ,  given  in  Equation  (3-42)  the 
best  derived  coefficients  are 


y*  =  1.3682 

6932 

2362 

r  =  -0.1268 

8142 

8338 

M  =  -0. 0024 

3108 

1471 

0^  =  -0. 1583 

3565 

5509 

Cq0  =  0.1014 

5485 

8689 

C^Q  =  0.0160 

1398 

7763 

Coi  =  -0.0904 

3332 

1481 

C^^  =  0.0574 

8768 

5173 

C20  = -0-0272 

7606 

9952 

Co2  = -0.0342 

1969 

6290 

Dq  =  0.0292 

5271 

6867 

=  0. 0059 

2537 

5028 
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APPENDIX  C  -  REFINEMENT  OF  THE  CARTESIAN 
TRANSFORMATION 


Section  3. 3  of  the  text  presented  a  numerically  efficient  and  compact  algorithm 
for  employing  the  Cartesian  transformation,  f  *  ,  defined  in  Section  2. 2  to  obtain 
the  data  base  coordinates  (x  ,  y)  of  any  datum  point  or  pair  of  datum  points 
on  the  (spherical)  Earth's  surface.  This  Appendix  discusses  the  rationale  for 
the  new  algorithms  and  provides  a  derivation  of  Equation  (3-42). 

C.l  OVERVIEW 

The  algorithm  presented  in  Reference  1  (Section  4,  Appendix  B)  already  pro¬ 
vides  the  essential  Information  needed  to  implement  f*  ,  but  this  study  has  found 
that  the  algorithm  may  be  made  more  efficient  computationally  and  better  suited 
to  computer  programming  by  the  following  modifications: 

•  Take  advantage  of  the  nondimensional  analysis  presented  in  Appen¬ 
dix  A  that  removes  the  need  to  compute  the  three-dimensional 
spatial  coordinates  of  datum  points  by  relying  on  simple  direction 
cosines  (1  ,  m  ,  n)  of  the  subsatellite  point  and  on  the  longii.u'L.es 
and  colatitudes  of  the  datum  points. 

•  Allow  for  the  Earth's  rotation  by  using  a  varying  node  (see  Appen¬ 
dix  A)  in  computing  the  direction  cosines  instead  of  performing  a 
separate  transformation. 

•  Reduce  by  about  half  the  computations  required  in  cases  of  nonzero 
scan  angles  by  viewing  the  scan  arc  as  a  central  angle  (see 
Appendix  A)  and  computing  the  data  base  coordinates  for  points  in 
symmetrically  placed  pairs  instead  of  treating  each  datmn  point 
independently. 

•  Remove  a  step  from  the  algorithm,  saving  computations  as  well 
as  yielding  desirable  numerical  properties  (see  Section  3.4),  by 
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introducing  a  change  of  variable  that  leads  to  a  new  formulation 
of  f*  .  The  results  of  this  step  are  summarized  in  Equations  (3-40), 
(3-41),  and  (3-42)  of  Section  3.3. 1.  The  mathematical  details  are 
presented  in  Section  C.  2  below. 

C.2  REFORMULATION  OF  f* 

Table  3-2  (Section  3. 3. 1)  establishes  the  conventions  that  define  variables  , 

T)'  ,  and  p  .  These  variables  also  appear  in  Table  4-1  of  Reference  1,  differing 
only  in  that  1  ,  m  ,  and  n  are  dimensionless  (see  Appendix  A),  whereas  a  , 
b  ,  and  c  (of  Reference  1)  arc  spatial  coordinates.  Mathematically,  4'  •  i 
and  p  behave  identically  here  and  In  Reference  1. 

la  Reference  1,  once  ,  r)'  ,  P  are  found  in  Table  4-1,  they  are  transformed 
to  (4 ,  n)  by 


«  P 

TJ  =  TJ'  P 


(C-1) 


where  P  =  r  /p  and  r  =  l/>/3"*  The  new  variables  £  and  ri  are  then  used 
o  o 

in  f*  to  obtain  the  data  base  coordinates  (x  ,  y)  by 


X  =  f*  iifV) 

y  -f*  (n,^) 


(C-2) 


The  transformation  given  by  Equations  (C-1)  can  be  eliminated  by  defining  new 
variables  a  and  j9  by 


“  =  r/p 

/3  =T?»/p 


(C-3) 
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Comparison  of  Equations  (C-1)  and  (C-3)  shows  that  the  new  variables  a  and 
eliminate  the  computation  of  the  intermediate  quantity  P  and  reduce  the  trans¬ 
formation  to  simple  ratios  of  previously  known  quantities.  The  introduction 
of  a  and  /3  will  now  be  completed  by  deriving  f’*'  (a ,  /3)  •  the  transformed 
equivalent  of  f*  (4  >  t?)  •  which  was  given  in  Section  2. 2  and  Reference  1  as 


f*  =  y  *4  +  (r^^  -  6*  +  6*  (r^^  - 


I  2  2\  V  *  ^2i  2j‘ 

j»0 


where  r  =  lA/s”  and 
0 


y*  =  i 

T  y 


6  = 


2r 


H*  = 


^  ^  1  ^3  -  2y*  -  p*  -  2r^  6*) 


<y  t  =  I  II 

y+  r  6 
o 


p+  2r^6 
o 


2r 


1  2 
r 

0 


(C-5) 


(Vj*  -Y*) 


_  6* 
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and  y  ,  n  t  6  ,  it} ,  and  the  and  d*  have  the  numerical  values  given  in 
Section  2. 2. 

First,  Equation  (C-4)  is  expressed  in  terms  of  4'  »  •  and  p  by  using  Equa¬ 

tions  (C-1) 


Ml'.  r}\  p)  =  y*r^  €'/p  +  r^  (4»/p)^ 


+  r^  (4»  Tj'^V)  [r^  -  r^  (C'/p)^]  U  +  6^  [r^  -  r^  (|'/p)^] 


j>0 


(C-6) 


+  r 


(|’/p)^[r^  -  r^  (|»/p)^] 


+[^0  -  (I'/P)^]  (C’/P)^  , 

'■  •'  iiO 


Neifit,  a  and  /3  are  introduced  as  defined  in  Equations  (C-3),  and  powers  of 
r^  are  canceled  in  Equation  (C-6),  where  appropriate,  to  obtain 


f*(o,  /3)  =  r^y*a  +  r^  (1  -  y*)  o‘ 


+  a 


0^(1-  a^)  r® 


6*  +  6*  r^  (1  -  a^) 
1  o  '  ' 


^  2  *  2i+2J  2i  -2j 

iiO 

+  (1  -  a^)  r®  j^w*  +  r^  (1  -  d*  r^‘  a^j 


(C-7) 
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Equation  (C-7)  may  be  simplified  ultimately  by  introducing  the  definitions 
(Equations  (C-5))  of  6*  ,  co*  ,  and  6*  to  obtain: 


f*(a,  ^)  =  r  +  r  (1  -  y*)  (1  -  a^)  r® 


,4  '^zl  4  ■  „4)*^‘“ 

2r  r  \  r  2r  / 

o  o  \  o  o  / 


iaO  ^ 

jso 


(C-8) 


^  3  2  6 

+  a  (1  -  a  )  r^' 


,3-  2y*  -M*  -  2r^  [oi*  -  H*)/2t^ 


2r 


.  2,,  2,  Vj*  2i  2ii 

+  r  (1  -  a  )2^d  r  « 

UO 


Equation  (C-8)  is  simplified  by  factoring  the  expression  by  r  and  carrying 
r  througli  the  terms  in  box  brackets  to  obtain 


f*(a,  8)  =  <y*a  +  (i  -  y*)  (i  -  a^)  j(p*  -  p*) 


+  (yj  -  y  -  YOil  -  M*))  (1 "  +  <1  -  18^)  23  % 

JiO 

+  (1  -  ^(3  -  2y*  -  +  >*♦  -  Mj)  +  (1  “ 

UO 


(C-9) 
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Three  new  constants,  M  ,  F  ,  and  ,  are  introduced  to  replace  6*  ,  on*  ^ 
and  through  the  following  definitions 

M  = 

T  =  y*-y*  (C-10) 


By  comparing  Equation  (C-10)  with  Equation  (C-5),  it  is  seen  that  these  new 
constants  are  considerably  simpler  in  form. 

Noting  that  r  is  strictly  a  constant,  r  =  l/\/T ,  new  coefficients  C,,  and  D, 
o  o  ij  i 

are  defined  as  follows 


j.2(i+j+3) 

o 


D 


i 


r2(i+3) 

o 


(C-11) 


Using  Equations  (C-10)  and  (C-11)  in  (C-9),  the  following  is  obtained 


f*(a,  /3)  =  r  |y*a  +  (1  -  y*)  (l  -  a^) 


M 


+  (T  -  M)  (1  -  a^)  +  (1  -  ^^)  Yi  C. . 

iiO  ^ 
j^o 


3  2,  fr.  ^  2,  „  2il  I 

a  (1  -  a  )  O.  +  (1  -  a  )  2^  D  o 

L  i^O  "  Jl 


(C-12) 
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This  completes  the  derivation  of  the  transformed  expression  for  f*  .  Equa¬ 
tion  (C-12)  is  given  as  Equation  (3-42)  in  Section  3.3. 1  and  its  numerical  ad¬ 
vantages  are  discussed  in  Section  3. 5. 


APPENDIX  D  -  APPROXIMATION  OF  THE  DIRECT  AND 
INVERSE  TRANSFORMATIONS 


This  zqipendix  presents  tables  of  coefRcients  for  the  approximate  forms  of  the 
transformations  f  and  f*  as  presented  in  Section  3.4.  These  tables  are  pre¬ 
sented  for  one  octant  of  one  face.  Each  table  contains  977  sets  of  coefficients 
(ordered  by  column),  determined  by  using  intersection  points  of  a  32  x  32  tri¬ 
angular  grid.  Each  set  of  coefficients  represents  the  ^proximate  form  of  the 
transformation  within  one  element  of  the  32  x  32  grid.  Each  of  these  elements 
corresponds  to  one  64  x  64  element  data  face  record  at  the  1. 2  nautical  mile 
scale. 

Figure  3-13  illustrates  the  table  octant  (octant  1)  and  the  seven  other  octants. 
Table  D-1  is  used  to  transform  from  (x  ,  y)  to  (^  ,  ri)  and  Table  D-2  is  used 
to  transform  from  ({•?))  to  (x  ,  y)  .  The  steps  outlined  below  are  followed 
to  transform  from  (x  ,  y)  to  (^ ,  ?)): 

1.  Add  2048.0  to  X  if  X  is  less  than  2048.0 

2.  Add  2048.0  to  y  if  y  is  less  than  2048.0 

3.  If  y  is  greater  than  x  ,  interchange  x  and  y  values 

4.  Truncate  the  x  and  y  values  obtained,  divide  each  by  32,  and  add 
one,  use  these  integers  and  IX  and  lY  indices  and  look  up  the  co¬ 
efficients  and  number  of  terms  in  the  table 

« 

5.  Evaluate  f  (x  ,  y)  to  obtain  $ 

6.  Reverse  steps  3,  2,  and  1  to  place  (  in  its  proper  octant 

7.  Interchange  x  and  y  and  repeat  steps  2  through  7  to  obtain  r) 

Alternatively,  steps  1,  2,  3,  and  7  may  be  avoided  by  constructing  a  iull 
64  X  64  table  containing  redimdant  data. 
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In  Tables  D-1  and  D-2,  IX  and  lY  (IP  ,  IE)  are  the  x  and  y  (^  ,  t))  indices  of 
the  table,  and  are  the  coordinates  of  the  origin  of  the  area  of 

definition  for  the  coefficient  set,  NTERMS  is  the  nunaber  of  terms  to  be  eval¬ 
uated,  a  ...  a.  are  the  coefficients,  and  RESID  is  the  largest  residual  found 
o  0 

in  the  evaluation  of  the  coefficients  (in  elements)  as  illustrated  in  Figure  3-12. 


D-16 


APPENDIX  E  -  LISTINGS  OF  TEST  AND 
EVALUATION  PROGRAMS 


E.l  SUMMARY 

This  appendix  contains  listings  of  FORTRAN  programs  used  for  testing  and 
evaluating  the  transformations  in  this  report  and  plotter  programs  used  for 


generating  the  world  maps . 

_ Listing _  Figure 

Program  used  for  evaluating  the  coefficients  of  f*({,t?)  E-1 

Program  used  for  verifying  the  f*  transformations  E-2 

Program  used  for  verifying  the  polar  transformation  E-3 

Program  used  for  evaluating  coefficients  of  at^roximation 

polynomials  E-4 

Program  used  for  generating  world  maps  E-5 

Alternative  subroutines  for  E-5  using  the  polar  transformation  E-6 
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